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Abstract 

After reviewing the basics of the cavity method in classical systems, we show how its quantum 
version, with some appropriate approximation scheme, can be used to study a system of spins with 
random ferromagnetic interactions and a random transverse field. The quantum cavity equations 
describing the ferromagnetic-paramagnetic phase transition can be transformed into the well-known 
problem of a classical directed polymer in a random medium. The glass transition of this polymer 
problem translates into the existence of a 'Griffiths phase' close to the quantum phase transition of 
the quantum spin problem, where the physics is dominated by rare events. The physical behaviour 
of random transverse field ferromagnets on the Bethe lattice is found to be very similar to the one 
found in finite dimensional systems, and the quantum cavity method gets back the known exact 
results of the one-dimensional problem. 
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I. INTRODUCTION 



The cavity method has been developed to study classical frustrated spin systems: spin 
glasses. In recent years its range of application has broadened a lot as it was applied with 
some success to hard computer science problems and also to some quantum problems. This 
paper aims at giving some background on the cavity method, both classical and quantum, 
and explaining how it can be used to study a quantum problem, a ferromagnetic spin system 
in a random transverse field, following the recent work of P'i]. This is a problem which has 
been studied a lot in one dimension; we analyze it here using the cavity method which is a 
mean-field type method better suited for large dimensional systems. However we find that 
the physics is rather similar to the one found in the one-dimensional case, in the sense that 
rare events (rare sites with anomalously small transverse fields) play a major role in the 
neighborhood of the quantum critical point. This leads to very large spatial fluctuations of 
the spontaneous magnetization in the ferromagnetic phase, and a broad distribution of the 
local magnetic susceptibilities in the paramagnetic phase. These effects, which are typical 
of the 'Griffiths phase' described in low dimensions, are found here through the use of an 
auxiliary problem of directed polymers in random media: the rare-event dominated regime 
of the ferromagnetic-paramagnetic phase transition is related to the glass transition of the 
polymer problem. 

This paper is organized as follows. Sec.|I| gives an introduction to the cavity method in 
classical spin systems. Sec.[TTl| defines the random transverse field ferromagnet and gives some 
background on this problem. The naive mean field approach to this problem is explained 
in Sec.|TV]. Sec. |V| introduces the quantum cavity method and discusses some systematic 
approximation schemes which can be used in this context. These approximation schemes 
are tested in Sec.|V| by applying them to a very simple, and very well understood, prob- 
lem: the pure ferromagnet in a transverse field. Sec.|Vl| describes the use of the quantum 
cavity method, with a 'cavity-mean-field' approximation, to study the phase diagram of a 
ferromagnet in a random transverse field. It derives the phase transition line in the plane 
disorder-temperature, and describes the main properties of the low-temperature ferromag- 
netic and paramagnetic phases, in particular the large fluctuations of the order parameter 
and of the susceptibility. This section ends by a study of the one- dimensional case, where 
the cavity method gets the known exact results, and a discussion of the higher- dimensional 
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cases, for which the Bethe lattice analysis provides a useful mean-field approximation. A 
short summary is given in Sec. |Vlli| . 



II. THE CLASSICAL CAVITY METHOD 
A. The Sherrington-Kirkpatrick model 

The cavity method has been introduced twenty five years ago in the context of spin 
glass mean field theory!. The solution of the fully connected Sherrington-Kirkpatrick (SK)S 




model had been found a few years before by G. ParisiH, who had proposed, within the replica 
approach, an inspired Ansatz of replica symmetry breaking (RSB). The cavity method gave 
an alternative solution, with identical physical content and results to the replica approach, 
but which was much more transparent. Instead of the mysterious replicas, it used a purely 
probabilistic approach, based on three assumptions which had been found to be hidden in 
the form of Parisi's RSB Ansatz: the existence of many pure statesi, ultrametricityi, and 
the exponential distribution of free energies of the pure states!. This approach has turned 
out to be very powerful and has served as a starting point of recent works by Talagrandi 
and Guerra0 which established rigorously the exactness of the free energy obtained by the 
RSB/cavity method (although the ultrametric property has still not be proven). 

In a nutshell, the cavity method as it was applied to the SK model consists in assum- 
ing some structure for a A-spin system s±, . . . , Sjv, and checking self consistently that this 
structure is reproduced when one adds a new spin s and goes to a N + 1 spin-system. 
Within the simplest replica symmetric (RS) approximation, the main assumption is that 
the correlation functions are small. The local field on sq then has a Gaussian distribution, 
with a width which must be reproduced self-consistently and gives the Edwards-Anderson 
order parameterS. This approach gives the RS solution, which is correct at temperatures 
above the de-Almeida Thouless (AT) line§. It breaks down in the spin glass phase, when 
the proliferation of pure states induces a non-trivial correlation between distinct spins. In 
order to study the spin glass phase, one must explicitly assume the existence of many pure 
states@0 with free energies distributed as a Poisson process with exponential distribution!. 
Within each state, the correlations are small and the cavity method can be applied. But 
the addition of the new spin s creates a free energy shift which depends on the state: the 
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free energies are reshuffled, and the selection of low free energy states (by the Boltzmann 
weight) biases the local field distribution. The whole cavity RSB approach relies on control- 
ling these crossings of the free energies of the states, and checking the self consistency of the 
basic hypotheses. The results are exactly identical to those of the replica approach. 



B. Finite connectivity spin glasses 

It was soon realized that it would be very instructive to go beyond the SK model by 
studying spin glass models in which each spin has only a finite number of neighbours. In 
order to introduce mean field models with this property, the usual approach is to study 
the Bethe lattice, which is usually defined as the interior of a Cayley tree: taking a large 
Cayley tree of depth M, one studies the interior part up to depth L. The Bethe lattice is 
then defined by the double limit Uhil-hx, liuiAf^oo (see Fig. |l]). This procedure is fine for 
ferromagnetic interactions, but it raises a problem when one studies spin glasses where a 
crucial ingredient is the frustration due to loops in which the product of exchange couplings 
is negative. A simple procedure is to induce frustration by fixing randomly (or imposing 
some type of local fields on) the spins at the boundary of the tree. However most of the 
physics is then put by hand through this boundary condition, which involves a finite fraction 
of the total number of spins. For instance fixing the boundary spins to ±1 independently 
with probability 1/ ■EE gives a system which does not have a genuine spin glass phase 
(technically it is is always replica symmetric). It has been found recently how the boundary 
conditions should be fixed in order to get RSB: the correct procedure involves the process 
of broadcast and reconstruction^ and creates subtly correlated boundary fields. 

An alternative construction is to define the Bethe lattice spin glass through a random 
regular graph, taken uniformly from all the graphs where each spin has K+l neighbours (see 
Fig. [I]). Large random regular graphs with K ^> 1 vertices have a locally tree-like structure 
like the interior of a Cayley tree, as the typical size of a loop is of order log N/ log K. These 
large loops create the frustration without having to play with boundary conditions. The 
study of spin glasses on random graph with a finite connectivity started with Viana and 
Bray@. The solution of spin glasses on random regular graphs within the RS approximation 
was found relatively early&il, but it took fifteen more years to understand how to handle 
RSB effects§, and this has been done so far only at the level of one, or at most two, steps of 
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FIG. 1. Left: the usual construction of the Bethe lattice. Right: a (small) random regular graph. 
The local structure in large random regular graph is tree-like, as in the center of a Bethe lattice. 



RSB. Here we shall briefly review how the cavity method can be used to study this problem, 
mostly at the RS level, as this gives the building blocks of the quantum cavity method which 
we use below. 

Let us study the equilibrium properties, at inverse temperature (3, of a general Ising spin 
glass problem with a random field on a Bethe lattice, denned by the Hamiltonian 

H = - ^ JijSiSj - 2J BiSi (1) 

<ij> i 

where < ij > are the links of a random-regular graph of degree K + l. The idea is to use the 
local tree-like property for rooted trees. Let us pick up one spin, say Sq, and delete the edge 
between and one of its neighbours. Up to any finite distance, Sq is now the root of a tree. 
Let us denote by ho the local 'cavity' magnetic field on the root sq due to this tree, so that 
(sq) = tanh(/3/io). Denote by s±, . . . , Sk its neighbours on the tree, and by h\, . . . , hk their 
cavity fields (so that for instance the magnetization of s\ in the absence of sq is tanh(/3/ii)). 
The Boltzmann measure on s is then 

K 



P(so) = ^5>xp(/3 



{si} 



B s + ^(JoiSoSi + hiSi 



i=i 



(2) 



By doing explicitely the sum over the spins Sj, one gets the recursion 

1 K 

h = B + — atanh(tanh(/3 J 0i ) tanh(/3/ij)) (3) 

which is illustrated pictorially in Fig.||]. One way to describe it is to build the tree rooted 
in s by considering the K trees rooted on si, . . . , Sk (in absence of s ) and merging them. 



FIG. 2. The recursion relation for the local magnetic field on the root of a rooted tree, here for the 
case K = 2. 

These recursion 'cavity' equations are always correct on a tree. On a random regular graph 
they are correct whenever the joint probability of the spins si,...,8k in the absence of 
So factorizes into a product exp (/3/ijSj) /(2cosh(/3/ij)). Two conditions are needed for 
this absence of correlations of neighbouring spins. One is the local tree-like property which 
guarantees that any two neighbouring spins Sj, Sj, with i, j G {1, . . . , K}, are far away on 
the graph obtained by eliminating the site sq and all edges connected to it. The second 
one is that the correlations decay at large distance, so that the correlations between two 
neighbours Sj and Sj (which are typically at distance logiV/ log-fT when sq is absent), vanish 
in the large N limit. This is true when the system is in its paramagnetic phase. It is not 
correct in a ferromagnetic phase or in a spin glass phase because long-range correlations 
develop. However, if one is able to restrict the measure to one single pure state (where by 
definition correlations decay at large distance), then the recursion relations are correct. 
These recursion relations can be used in two main ways: 

1- On a given sample, one can see these equations as a set of self-consistent equations 
relating cavity fields: there are two such cavity fields for each edge < ij > of the graph, one 
is the field on site i in the absence of j (case where i has been chosen as root), the other one 
is the field on site j in absence of i. All these fields are related by the cavity equations. In the 
large K limit these equations reduce to the 'TAP' equations^. These RS cavity equations 
are also well known in computer science under the name of 'belief propagation' equationsil. 

2- One can also study the statistics of cavity fields. Suppose that the local magnetic fields 
Bi are independent random variables drawn from a distribution P ext (B), and the exchange 
couplings are independent random variables drawn from a distribution p(J). At the RS 
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level, one expects a single solution for all cavity fields on a given instance. Then one can 
define the distribution of cavity fields P(h) when one picks up an oriented edge randomly. 
This is the natural order parameter for this problem. Its knowledge provides the full solution 
of the problem when the RS hypothesis is correct. Technically, finding P(h) is usually done 
by a 'population dynamics' method first introduced in [H], and developed in the present 
context by [0]. Its idea is the following. The update equation induces a self-consistency 
equation for P(h) which can be written as: 



K / -, K 



P(h) = J dBP ext (B) f[ J dhiP(hi)dJip(Ji)6 \ h - B - ^ atanh(tanh(/3 J t ) tanh(/3/ii)) 

(4 ) 

The population dynamics method represents P(h) by a population of M > 1 fields 
hi, ... , Hm- This population is updated by a Monte-Carlo-type process in which, at each 
iteration, one does the following operations: 

• Choose K indices 7*1, . . . r# G {1, . . . , M} randomly uniformly. 

• Generate K independent couplings Ji, . . . Jk from p(J), and a field B from P ext (B). 

• Compute h = B + 4 Yl^=i atanh(tanh(/3Jj) tanh(/3/i r J 

• Choose an index j G {1, . . . , M} randomly uniformly, and replace in the population 
the value hj by the new value h. 

The convergence of this method when the number of iterations is large enough can be 
checked by monitoring for instance moments of P(h) (the r-th moment is evaluated as 

enough the population will give a good approximation of 

P(h). 

Note that the RSB situation is much more complicated: on a given instance there are 
many solutions to the system of equations relating the cavity fields. On a given oriented 
edge (say looking at the cavity field on i in absence of j) the cavity field can take many 
values depending on the solution one takes. Doing the statistics over all solutions defines a 
distribution of fields on each oriented edge, and the order parameter is the distribution of 
these distributions when one picks an edge at random. In the following quantum problems 
we shall use only the RS version. 



III. RANDOM TRANSVERSE-FIELD FERROMAGNETS 



The problem of ferromagnets in random transverse fields has received a lot of attention 
in recent years. These systems provide relatively simple examples of disordered systems 
displaying a quantum phase transition^. They are described by the Hamiltonian: 

11 ' (5) 

* (tf) 

where a z , a x are Pauli matrices. The exchange couplings are independent ferromagnetic 
interactions drawn from a distribution p(J) which has support on J > 0, and the random 
transverse fields are independent variables drawn from a distribution 7r(£). The lattice 
structure is described by the set of pairs (ij) which appear in the above sum; we suppose 
that this lattice is homogeneous, and each spin has exactly z — K + 1 neighbours. For 
definiteness we shall use a model where the distribution of couplings and fields are given by: 

p{J) = *9{J)6{2/K-J) ■ tt(0 = \e{i)9{h - £) , (6) 

where 9 is the Heaviside function. The choice of the width of the J distribution fixes the 
energy scale such that the mean-field critical temperature in absence of disorder has T c — 1. 
The Hamiltonian Hp thus depends on a single parameter h which measures the degree of 
disorder. 

The one dimensional case, which is intimately related to the two-dimensional classical 
model of MacCoy and Wull, has been solved in great detail by FisheiH using the strong 
disorder Ma-Dasgupta-Hu decimation procedure!^!. This solution has emphasized the im- 
portance of Griffiths singularities, which manifest themselves most notably in the fact that 
the average susceptibility differs from the typical one. The method itself has a broad range 
of applications^. These results have been confirmed numericallyil'0, and numerical simu- 
lations of the same problem in dimension two, both through the strong disorder decimation 
and through Monte-Carlo simulations, also point to the same kind of physics as in one 
dimensionilUl. Recent results also find the same behaviour in dimension three and higherU 

Using the method introduced in |0i], we shall now study the Bethe approximation for 
the random transverse-field ferromagnet described by (H). We thus assume that the spins 
are on the vertices of a random regular graph of degree z — K + 1 and we want to study the 
phase diagram of this problem as function of the temperature T and the amount of disorder 
h. 



8 



IV. MEAN FIELD 



A. Naive mean field 

For large z one may be tempted to use the mean-field approach, in which Hp is replaced 
by Hmf = Tlii—tii&i ~ Baf) where B is determined self-consistently from the equation 
B = ^2j(<Jj). At temperature T = 1/(3, the self consistent equation for B obtained at 
z — > oo is: 

B = z j d^(0 J dJp(J) ^ / _^_ tanh(/3 v ^+~^) . (7) 

At finite temperature there is a transition between a large disorder paramagnetic phase 
where B = is the only solution of this equation, to a ferromagnetic phase characterized by 
a non-zero field, B > 0. The critical temperature T c is given by: 

if^y.i. ( 8) 

"> Jo s 

At zero temperature this naive mean field approach predicts that the system is always in its 
ferromagnetic phase. 

While this approach is correct when z = oo, its conclusions are qualitatively wrong for 
any finite z. 



B. Bethe-Peierls approximation 

The naive mean field prediction that the system is always ferromagnetic at zero T is 
wrong; the reason is that the mean field approach does not take into account properly the 
rare events. Rare sites j where the energy |^ | <C 1 can be easily polarized in the x direction, 
and they play a key role in the establishment of the ferromagnetic order. This effect, which 
has been first discussed qualitatively in leads to a zero temperature quantum phase 
transition to a paramagnetic state at a critical value of disorder h c . In order to study this 
quantum transition in details, we shall use the quantum version of the cavity method. 
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V. THE QUANTUM CAVITY METHOD AND ITS PROJECTIONS 



A. The RS quantum cavity method 

Let us study the transverse-field Ising Hamiltonian (|5|) on a Bethe lattice defined as a 
random regular graph. The partition function Z = Tr e _/3H/ can be expressed with the 
Suzuki- Trotter representation. Using M imaginary time steps, we introduce at each time 
step a decomposition on the eigenvectors of the operators crj. This describes the system by 
the classical time trajectory of each spin, <7j(t) G ±1 , where t = 1, . . . , M, and <7j(t) denotes 
the eigenvalue of cr? at time t. Then 



Z= lim V e~ pHsT (9) 

where: 

HsT = -Ji E E J**i(t)<ri(t) -EE r^V^ + 1) ( 10 ) 

t (ij) t i 

and Ti = ^logcoth 

This Hamiltonian acts on spin trajectories 5j = {(Tj(t), t = 1, . . . , M}. Each trajectory 
can be seen as a variable <7j taking 2 M possible values, and Hgr defines the interaction of 
these iV variables. The crucial point is that this Hamiltonian involves interactions between 
pairs of trajectories which are neighbours on the Bethe lattice. This is a locally tree-like 
graph and therefore one can use the cavity method to study it. The RS cavity method for 
this problem was introduced in |@|. It is obtained from the classical cavity method described 
in Sec. [IIB| by using spin trajectories instead of Ising spins. Take a branch of the Bethe lattice 
rooted on spin 0, and define the probability distribution of the time-trajectory of this spin 
as %pQ [<jq\. Define analogously, for each of the K spins i which are the first neighbors of 
on the rooted tree, their time trajectory as ipi [Bj\. Then one can write the mapping that 
generates, from {ipi [<?j]}, the new ipo [<?o]: 

ipQ [a ] = C ^ [°i] • • • $K [&k] 

<71,...,<TJf 

ex P \Jf E E J oMt)ai(t) + /? E T oVo(t)°o(t + 1) (11) 

\ t (i) t J 

The population dynamics method of Sect. HB| can be used to study this cavity recursion: one 
must represent each of the ipi by a sample of spin trajectories. This has been done for some 
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problems in , but this approach tends to be rather heavy numerically. In particular, in a 
disordered system like the random transverse-field ferromagnet, the natural order parameter 
is the distribution of the functions ipj when the site j is drawn randomly. This is a rather 
complicated object which must be represented as a population of populations. 



B. The projected cavity mapping 

While the exact RS cavity mapping is in principle doable, the difficulty in its numerical 
resolution may make it difficult to get a clear understanding of the physical results. For 
this reason one may want to develop approximate versions of it. One such approximation, 
the projected cavity mapping introduced in p\, consists in using a local distribution ipi [<7j] 
which takes the special form, parameterized by one single number B{. 




J2^(t) + f3T l J2^( t )^( t + 1 n ■ ( 12 ) 

t t J 



ipi [<Ji (£)] = Cexp 

Injecting this form into ( |TT|) we obtain a ipo which does not have the form (]12|) . However 



one can project it back on a i/jq in the right subspace described by (J12[) as follows. From the 
obtained i[)q one can compute (ctq) and deduce from it the effective parameter B such that 
the distribution C exp a (t) + /3T Y^ t o"o(t)cr (t + 1)) gives this value of 

This projected cavity mapping is actually more easily understood directly in terms of 
quantum Hamiltonians, without going to the Suzuki- Trotter formalism. One studies the 
properties of a spin in the rooted graph where one of its neighbors has been deleted, 
assuming that the K remaining neighbors are uncorrelated. Each of these neighbours, when 
it is at the root of the subtree obtained by deleting the link to spin 0, is supposed to be 
described by a Hamiltonian 



Ht av = -to - B i0 f (13) 
The system of spin and its K neighbors is thus described by the local Hamiltonian 

K 

H = -totf - (&*f + B i°t + J<*°o<>i) ( 14 ) 

i=l 

By diagonalizing this 2 K+l x 2 K+1 Hamiltonian, one can compute the induced magnetization 
of spin 0, mo = In order for this approach to be self-consistent, this magnetization 
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should be equal to the one obtained from the Hamiltonian — ^af — Bo<jf. This means that 
Bq is obtained by solving the equation: 



B n 



tanhpJ$ + B$ = mo (15) 



The projected cavity mapping thus uses a single number Bi to describe the trajectory dis- 
tribution ipi. It leads to a mapping which gives the new cavity field Bq in terms of the K 
fields Bi on the neighboring spins. This mapping induces a self-consistent equation for the 
distribution of the B fields, which is the natural order parameter in this context: the para- 
magnetic phase has P(B) = 6(B), and the ferromagnetic phase is described by a non-trivial 
P{B). The population dynamics method can be used to find this distribution. 



C. The cavity-mean-field approximation 

The projected cavity mapping is numerically much simpler than the full RS cavity method 
using spin trajectories. It also allows to address a broader range of questions, like those 
related to the real time dynamics and the spin relaxation, which are not easily accessible 
within the Suzuki- Trotter formalism!'!. Still its practical use is a bit slow numerically as 
one must diagonalize the K + 1 spin Hamiltonian flU] ) each time one wants to generate a 
new field Bq in the population. 

It has turned out usefulS to use one more step of approximation in order to obtain an 
explicit mapping, similar to the one found in the classical problem (^). This can be done 
using a mean field approximation in order to compute the magnetization m from the cavity 
Hamiltonian ([14]). In this ' cavity- mean-field approximation ', one approximates the cavity 
Hamiltonian acting on spin by 

K 

K V ' MF = - < E J M) (16) 

i=l 

This implies that B = 5^ i=1 Joi(crf), giving the explicit recursion equation relating the B 
fields: 

K B l 

i=i v B i +£i 

This recursion induces a self-consistent equation for the distribution P(B) which can now 
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be written explicitly: 



K / K 



P(B) = I n md^dJ 0l P(B t )p(J M^)] Si B - "7=|=f + $ 

J i=i \ i=i V B i +Qi 



(18) 



VI. A SIMPLE TEST OF THE PROJECTED CAVITY RECURSION: THE PURE 
FERROMAGNET 

It is useful to test the validity of the two levels of approximation described above (the 
projected cavity mapping and the cavity-mean-field approximation) versus precise results 
using the full numerical sampling of Suzuki- Trotter trajectories. This can be done in the 
case of a pure ferromagnet in a uniform transverse field, defined by p( J) = 5{J — 1) and 
7r(£) = 5(£ — h). We thus study the problem defined by the Hamiltonian 

H PF = -hj2^-J2 a ^ ( 19 ) 

on a random regular graph where each site has z = K + 1 neighbours. The exact cavity 
study in terms of continuous-time spin trajectories has been done in the case K = 2 by [H], 
we shall use their results as a benchmark in order to test the accuracy of the three levels of 
appr oximat ions : 

1: Naive mean field: The Eq.(^) for the spontaneous field B in the x direction becomes: 

B 



5 = (K + l)- ? =_tanh(/3v / ^ + ^) . (20) 

2: Cavity-mean-field: The recursion relation (|T7D now becomes a self-consistent equa- 
tion for the spontaneous field B in the x direction (which is site-independent): 

B 



B = K-j=== tanh(/37)^T^), (21) 

it differs from the naive mean field only by the replacement K + 1 — > K. 

3: Projected cavity mapping: The system of spin and its K neighbors is described 
by the local Hamiltonian 

K 

Ho = -hoi - E ( ha i + Bo * + °o°i) ( 22 ) 



i=l 
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The self-consistent equation for the spontaneous field B in the x direction is given by 

Tr(a$e-P Ho ) B 



m = A o L = tanh(3Vh 2 + B 2 . (23) 

Tre~^ ^h 2 + B 2 V ' 

These traces involve matrices of size 2 K+1 , they are easily computed numerically when K is 
not too large. 

The phase diagram obtained with these three approximations is plotted in Fig.|3], and 
compared to the exact one of |H|. The projected cavity mapping gives a rather precise 
approximation of the phase diagram, much better than the one obtained by naive mean 
field or by the well known static approximation (see 

pi]). 

A slightly better result can be 
obtained by using the parameterized form (12) variationally inside the Bethe free energy^. 

One may also use the cavity method as a mean field approximation to study the finite 
dimensional problem. It turns out that it is able to locate the transition point rather 
accurately. For a three dimensional case, the simulations through a cluster Monte-Carlo 
method of the Suzuki- Trotter formulation^ gives a zero-temperature transition occuring at 
h c = 5.16. The projected cavity mapping (with K — 5) gives h c = 5.28, and the cavity- 
mean-field gives h c = 5.00. In two dimensions, the numerical result of [H0| is h c = 3.04, 
the projected cavity mapping (with K = 3) gives h c = 3.22, and the cavity-mean-field gives 
h c = 3.00. 



VII. PHASE DIAGRAM OF THE RANDOM TRANSVERSE-FIELD FERRO- 
MAGNET 

Let us apply the cavity-mean-field approximation to the random transverse-field ferro- 
magnet. The whole problem reduces to finding the distribution of fields P(B) which solves 
the self-consistent equation flT8|) . It turns out that the solution of this equation is more 
subtle than the one of seemingly similar equations obtained in the classical case. 

In order to understand the properties of this solution, it is useful to go back to the cavity 
mapping fll7|). Let us iterate this mapping L ^> 1 times on the Bethe lattice. For L finite 
and iV — y oo the corresponding graph is just a rooted tree with branching factor K at each 
node and depth L. The field Bq at the root is a function of the K L fields on the boundary. 
In order to see whether the system develops a spontaneous ferromagnetic order, we study 
the value of B in linear response to infinitesimal fields Bi = B 1 on the boundary spins. 
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FIG. 3. Top: The phase diagram of the ferromagnet in a transverse field, on a Bethe lattice 
with K + 1 = 3 neighbours per spin. The plot shows the critical temperature versus transverse 
field, separating the low-T ferromagnetic region from the high-T paramagnetic one. The full red 



green long-dashed curve is the naive mean field result. The blue dashed line is the result from the 
cavity-mean-field approximation. The dotted purple curve is the result obtained from the projected 
cavity mapping. Bottom: Same plots with K + 1 = 6 (without the exact result), showing that 
the approximation going from the projected cavity mapping to the cavity- mean-field improves at 
larger K, as expected, and becomes rather accurate already for this moderate value of K. 

This is given by 



where the sum is over all paths going from the root to the boundary, and the product Ylk 
is over all edges along the path P. 



curve is 



the exact result obtained by continuous time spin trajectory population dynamicsEl The 





(24) 
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A. Phase transition and directed polymers 



The response H is nothing but the partition function for a directed polymer in a random 
medium (DPRM) on a tree, where there exists on each edge k of the tree a random energy 
Ek given by e~ Ek = Jfc(tanh(/3£fc)/£fc) and the temperature of the polymer has been set equal 
to one. The random energies are independent identically distributed random variables. We 
need to compute the large L behavior of S. As S is a partition function, this is naturally 
characterized by the free energy 



$ = lim (1/L)lni 

L—>oo 



(25) 



If $ < the effect of a small boundary field decays with distance: the spin system is 
paramagnetic. If $ > it is ferromagnetic. 

The general method for computing $ in DPRM on trees has been developed by Derrida 
and Spohn0. Their solution can be expressed in terms of the convex function 

1 



/(*) 



In 



x 



(26) 



Let us denote by x = m the value of x G [0, 1] where this function reaches its minimum.The 
free energy of the DP is then $ = f(m). The phase transition line in the (h,T) plane 
separating the ferromagnetic from the paramagnetic phase is thus obtained by solving the 
equation f{m) = 0. 

In the case of the random transverse-field ferromagnet one obtains 

'K l ~ xl 



/(*) 



In 



x 



1 + x 



, 2 1 , 
+ log - + - lo§ 

1 X 



T f h/T , Aanh(w) 
du I 



h 



u 



(27) 



from which one deduces the following phase diagram (see Fig.| 



When the field h is smaller than a value Hr the phase transition line is given by 
/(l) = 0. This gives back the naive mean field result of (|8|): h = j^ Tc duta,nh(u) /u. 
The value Jir, and the corresponding value Tr = T c {Hr), are obtained by solving the 
pair of equations /(l) = 0, /'(l) = 0. 

When h > h^, the phase transition is given by solving the equations f(m) = 
0, f'{m) = 0. In this regime the solution is at m < 1. This low temperature, 
large disorder part of the phase transition line departs from the naive mean field 
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prediction. In particular, at variance with the naive mean field result, one finds that 
there exists a zero temperature quantum phase transition at a critical disorder h c . 
This is obtained using the zero-temperature limit of /: 

f (x) = 1^ InK - 1 log(l - x 2 ) - log 0£\ , (28) 

and solving the two equations /o( m ) = fbi m ) = 0- Fig-|§ shows h c as function of K. 
In the large K limit fa c diverges, giving back the naive mean field result. However this 
divergence is very slow: one finds that m ~ 1 — 1/ \ogK and h c ~ elogK. 

The phase transition line has two parts: a small-disorder part on the left of the point 
R with coordinates (/ir, Tr), and a large-disorder part on the right of this point. The 
properties of these two regimes of the transition are very different. In order to understand 
this difference, it is useful to give a closer look at the DPRM problem using the replica 
method, as was first done in [0]. The DP partition function H, defined in (|24|), depends on 
the random quenched variables J n ,£ n . To compute the value of logS for a typical sample 



one needs to do a quenched average of logS over these random variables, denoted by logS. 
As a first step one might try to approximate it by the annealed average logH. This gives 
precisely the naive-mean field result. To go beyond this 'annealed approximation' one can 



introduce replicas, using InS = lim n _ s>0 (" ?1 — !)• The average of S n is obtained by a sum 
over n paths, 



E n 

P,,...,P„a=l 



(29) 



Lk&Pa 

One can then proceed by testing various hypotheses concerning the structure of the paths 
which dominate the replicated partition function (|29|) : 



• In the so-called replica symmetric (RS) hypothesis one assumes that the leading con- 
tribution to ( ]2"9"| ) comes from non-overlapping paths. This gives 

^ = exp(Ln/(l)) = (H)" . (30) 

The RS hypothesis is equivalent to the annealed approximation in this case. 

• The so-called one step replica-symmetry-breaking (1RSB) hypothesis assumes that 
the leading contribution to (|29|) comes from patterns of n paths which consist of n/m 
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FIG. 4. Phase diagram of the quantum ferromagnet in a random transverse field, as found with 
the RS cavity method, in the plane of the disorder strength h and the temperature T. The RS 
transition line (solid) is determined by the condition on the convex function: f(x = l,T,h) = 0, 
and separates the ferromagnetic (FM) state (f(x = l,T,h) < 0) from the paramagnetic (PM) 
state {f{x = 1,T, h) > 0). The point R on the RS line corresponds to the minimum of the convex 
function: f' x (x = l,T,h) = 0. It is the beginning of the RSB transition line, determined by the 
minimum of the convex function in the RSB (x* < 1) phase: f(x*,T,h) = and f^(x* ,T,h) = 0. 
The inset shows details of the RSB lines in the cases K = 3 and K = 5. The critical values of 
disorder h c (K), at which the zero-temperature quantum phase transition to the paramagnetic phase 
happens, for the two connectivity numbers K are: h c (K = 3) = 8.36 and h c {K = 5) = 10.28. 

groups of m identical paths, where the various groups go through distinct edges. This 
gives: 

W l = exp(Lnf(m)) (31) 



where f(x) is the function introduced in (p7|) . In the replica limit n — > 0, the parameter 
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FIG. 5. The critical value of disorder h c (K), as a function of the connectivity number K. 



m should belong to the interval [0,1]. It turns out that one should minimize the 
function f(m) over m G [0, 1] (the fact that one should minimize /, and not maximize 
it as one would have thought naively, is a well-known aspect of the replica method^!) . 
The low-temperature large-disorder part of the phase transition line corresponds to 
this 1RSB solution. 



B. Replica symmetry broken phase of the polymer and Griffiths phase of the 
ferromagnet 

This replica solution of our DPRM, with a RS phase at small disorder and a 1RSB phase 
at large disorder, fully agrees with the original travelling wave solution of [0]. So the point 
R on the paramagnetic-ferromagnetic transition line of the quantum random transverse-field 
ferromagnet corresponds to a 1RSB transition in the auxiliary DPRM problem. This is a 
glass transition of the polymer, similar to the one found in the random energy modeS0, 
where the measure on the paths condenses on a finite number of paths. A useful order 
parameter to describe this transition is 

" p keP 

This is nothing but the sum ^ p wp where Wp is the relative weight of path P in the measure 
(^|). Its inverse gives a measure of the effective number of paths contributing to S. In the 
RS phase one finds that Y vanishes in a typical sample. This can be understood from the 
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J k tanh(/3£ fc ) 



(32) 



fact that the DPRM partition function gets a contribution from a large number of paths 
(diverging in the large L limit). In fact the DPRM problem has a finite entropy in this RS 
phase. On the contrary, in the 1RSB phase one finds that Y is finite and fluctuates strongly 
from sample to sample. For instance, its first moments are given by p\: 

Y=l-m ; T i = ~(l-m) + ~(l-m) 2 , (33) 

3 3 

and the entropy density of the polymer vanishes in this whole 1RSB phase. 

The DPRM problem is an auxiliary construction used to study the phase transition in 
the original problem of the random transverse field ferromagnet. So the natural question 
is: what are the consequences of the 1RSB transition in the DPRM concerning the original 
quantum spin system? Let us go back to ( plf) , which connects the two problems. We consider 
a large tree of length L ^> 1. The condensation of the paths in the DPRM means that, if 
a small field is applied to the boundary, the value of the effective field Bq on the root is 
dominated by a finite number of paths to the boundary. So the root spin feels the effect of 
only a finite number of the boundary spins. The physics is dominated by rare events: rare 
boundary spins have a dominant influence of the central spin, and it is building up on these 
rare events that the system can or cannot develop a long range ferromagnetic order. This 
is very similar to the Griffiths phase, found originally in the one-dimensional case@, and 
argued to be present at least in low- dimensional systems. It is interesting to see that we find 
here a similar rare-event dominated phase, using the mean-field type Bethe approximation 
which is valid in the opposite limit of large connectivity, K 1. The next subsections 
will give a more detailed discussion of the connections between the results obtained with 
the cavity method and those known, or conjectured, for finite-dimensional problems. Let us 
first summarize the main observables that are accessible using the cavity method. 

One can study the practical consequences effect of the rare events on both sides of the 
phase transition, in the neighborhood of the phase transition line, on the right of point R. The 
study of |i] has discussed these effects in the context of the superconductor-insulator transi- 
tion. It is easy to translate them to the present situation of the ferromagnetic-paramagnetic 
transition. 

• In the ferromagnetic phase the order parameter (here we use as order parameter the 
typical local magnetic field B typ in the z direction, obtained as B typ = exp(logZ)) has 
an anomalous scaling and distribution. Imagine one enters the ferromagnetic phase 
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at a given temperature T by varying the disorder h. At large enough temperature 
the scaling is B typ ~ \Jh c — h, as one expects in a mean field system. Then lowering 
the temperature there is a region above Tr where the scaling is B typ ~ (h c — h)" 1 ^ . 
The exponent a(T) diverges when T — >■ Tr. In the 1RSB regime, i.e. for T < Tr, the 
typical field has an essential singularity in e G ^ hc ~ h \ In this regime the distribution 
of the local order parameter has a power law tail B^ p / B l+m . As m < 1, the average 
order parameter is divergent. 

• In the paramagnetic phase the local susceptibility fluctuates strongly from site to site 
and nonlinear effects are crucial. Consider the system in a very small uniform magnetic 
field B ext and look at the local order parameters Bi induced by this external field. The 
typical value of Bi scales linearly with B ext , but the moments B x with x > m, and 
in particular the mean B, are divergent at the level of the linear response to h: they 
behave non-linear ly, as B x ~ CB™ t . The distribution of the local susceptibilities 
Xi = Bi/B ext has a power law tail P(x) ~ C/x 1+m - 



C. The one-dimensional case 



The one-dimensional RTFF is the best understood case. The strong disorder decimation 
procedure of Ma-Dasgupta-Hu, developed by FisherEl, becomes exact at the transition point 
because of the existence of an infinite disorder fixed point. It can also be checked versus 
alternative methods like the mapping of the problem to free fermionsH'll, and the exact 
solution of McCoy and Wu@. A large corpus of results has been obtained by these methods, 
concerning the scaling behaviour and the Griffiths phase. The strong disorder decimation 
has turned out to be useful in a broad range of problems^. 

Let us mention only some aspects of the physics of the one-dimensional case for which 
one can draw a comparison between the exact results and those of the cavity method. 
The exact analysis of the problem has shown that the zero-temperature critical point is 
located at the width of the random transverse field h = h c such that log £ = log J, or 
explicitely: j <i£7r(£)log£ = J dJp(J)\ogJ. In our case this gives h c = 2. At the critical 
point the system displays activated dynamical scaling: the scaling exponent z c , relating the 
characteristic time scale r to the characteristic length scale £ through r = C£ 2c , is formally 
infinite. Instead the scaling is logr = Ct^. The Griffiths phase can be studied close to 
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the critical point. Defining 5 = (logh — log J) /(var [log h] + var[logJ]), where the overline 
denotes the average and var denotes the variance, this exponent has been shown to behave 
as z(5) ~ 2/5 + C, where C is a non-universal constant. This exponent z(5) also appears 
in the singular response of the system to an external magnetic field H: in the paramagnetic 
phase, the linear susceptibility diverges, and the magnetization behaves as M ~ \H\^ zt ^. 

Let us now present the results of the cavity method described in the previous subsection. 
The one-dimensional case is obtained by taking K — 1. Actually it is useful to keep K = 1+e 
at intermediate steps of the computation, and take the e — > limit in the end. At zero 
temperature, the minimum of the function f(x) is obtained at x = m = -y/e, and the critical 
value of h, obtained by solving f(m) = 0, is h c = 2(1 + l^fe). When e — > 0, h c goes to 
2 which is the exact result, and m — > 0. The exponent m is the one that controls the 
decay of the local field distribution (which goes as 5~( 1+m )) ; and the distribution of local 
susceptibilities. From the analysis of p], this exponent m controls the response to a uniform 
external field, M ~ |if | m . So one should identify m = l/z(5). At the critical point we have 
found m — > 0. This agrees with z c = oo and the activated scaling behaviour. Let us now 
study the behaviour close to the critical point, at h = h c (l + y). At K = 1 (or e = 0), one 
obtains m = y + 0(y 2 ). So the prediction of the cavity method is z = 1/y. It turns out that, 
when expressed in terms of the distance to the critical point 5, this gives the exact result 
z = 1/(25). 

The quantum cavity method is exact in one dimension, if one uses the full mapping of 



spin trajectories given in (11). It must thus reproduce the exact results. What we have 
seen is that even the simpler cavity-mean-field approximation gives a fair description of the 
problem, including the exact location of the critical point, the activated dynamical scaling, 
and the correct value of z close to the critical point. Actually it gives the exact correlation 
even for a finite length problem, as we now show. 
Consider a chain described by the Hamiltonian 

L L-1 

Hid = ~ 22 ~ Jn(J nK+l ■ ( 34 ) 



n=l n=l 



In order to see the onset of long-range order, one can fix that the end-spin is in the eigen- 
state of <t£ with eigenvalue +1, fix £^ = 0, and compute the so-called 'surface' magnetization 
at the other end-point in the ground state m\ = (0\af |0). By using the Jordan- Wigner trans- 
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formation of this problem to free fermions, Igloi and Rieger have shown thatE2l: 

-1/2 



mi 



L-l / n 



i+E n 

n=l \r=l 




(35) 



This formula, which is exact for a finite system, can then be used to study the critical 
behavior of the problem. For instance, [@] shows how to deduce from this formula several 
interesting physical properties: the typical and the average surface magnetization differ, the 
average critical magnetization scales like L _1//2 , the exponent (3 S giving the surface magneti- 
zation close to the critical point is equal to j3 s = 1, the correlation length exponent exponent 
v is equal to v = 2. 

Let us now use the cavity-mean-field method on this one-dimensional problem. At zero 
temperature, formula flT7| ) gives the recursion 

B„-i = J n _i^===, (36) 



which should be initialized with Bl = oo, £l = 0. It is easy to see that, iterating the 



mapping fl3"5D, and computing rrii = B x j \/ B\ + C,f, one obtains the exact formula ([35]). 
Therefore the cavity method, even using the simplified, mean-field mapping, gives the exact 
value of the surface magnetization, also for a finite system. Its predictions thus agree with 
the known results, including the values of the exponents f3 s and v. 

D. Dimensions two and above 

There is no exact solution of the two-dimensional RTFF problem. This problem has been 
studied mostly using the strong-disorder decimation procedureSilS, and also using cluster 
Monte-Carlo simulations of the 2+1 dimensional problem obtained in the Suzuki- Trotter 
representationillii. These works have mostly used a different normalization, in which p(J) 
is the uniform distribution on J 6 [0, 1], and 7r(£) is uniform on [0, h]. In this subsection we 
shall thus adopt this other normalization. The results obtained both by decimation and by 
Monte-Carlo indicate that the situation is similar to the one-dimensional case. In particular, 
the exponent z c is found to diverge, and l/z(6) goes to zero as 5 — > 0. In order to obtain an 
analytic approximation to this two-dimensional problem, one can repeat the cavity study 
on the Bethe lattice with K = 3. 



23 



Within the cavity- mean- field approximation, the function f(x) is found to be 



f(x) = — log 3 log(l — x 2 ) — log h. 



(37) 



x x 



It is minimum at x = m = .679, and the critical value of the field, obtained from f(m) = 0, is 
found to be h c = 12.53. This result disagrees with the strong-disorder decimation, and with 
the Monte-Carlo, in two aspects: i) the value of h c (with the Monte-Carlo it is foundiilll 



the fact that the value of m found in the cavity method does not go to zero at the transition. 

One may wonder to what extent these results depend on the approximation. If instead 
of doing the cavity-mean-field approximation one uses the more accurate projected cavity 
mapping, the value of h c is changed to h c ~ 7.5, but the value of m is unchanged. This is 
seen from Fig|| which shows the inverse of the typical value of the local field B, measured 
as exp(— In I?), plotted versus (h c /h) m — 1. The prediction of P] is that the inverse of the 
typical value of B should go to zero linearly in (h c /h) m — 1. Both data sets, obtained 
from the cavity-mean-field approximation and from the projected cavity mapping, find this 
behaviour, but the values of h c in both data sets are distinct (respectively h c = 12.53, as 
predicted from our analytic study with the DPRM, and h c = 7.5, found from the fit). 

The conclusion is that, while the non-universal critical value h c depends quite a lot of 
the approximation used in the cavity method, the value of m does not, and the statement 
that m ^ at the transition seems to be a solid result for the problem on the Bethe 
lattice with any connectivity K + 1 > 3. Notice that this result still implies that there is a 
strong hierarchy of local magnetic fields: the average field diverges, and so does the average 
susceptibility. But it is not an infinitely strong hierarchy as happens when m = 0. It should 
be noticed that our best cavity analysis, using projected cavity mapping, still involves some 
approximation, and it would be interesting to carry out the full Suzuki- Trotter trajectory 
mapping in order to get the exact solution for the Bethe lattice. At present, our approximate 
treatment gives m > at the transition for the Bethe lattice problem a result which disagrees 
with a recent study from strong disorder decimatiorfi It should be noticed that our cavity 
approximations shoud become exact in the large K limit of the Bethe lattice. Assuimng that 
the result z c ^ oo is correct for the Bethe lattice, this raises the question of the existence 
and value of a critical dimension above which the finite dimensional problems have z c ^ oo. 

It would thus be interesting to revisit the numerical results in d = 2 and d > 2, and on 



around h c ~ 4.2, with the decimation it is 

foundOT 

around h r ~ 5. 




3), ii) most importantly, 
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FIG. 6. Phase transition in the RTFF model with K = 3, a uniform distribution of couplings J on 
[0,1], and a uniform distribution of fields £ on [0, h]. The typical value of the local field B, measured 
as exp(— lni?), is plotted versus (h c /h) m — 1. The crosses are obtained with the cavity-mean-field 
approximation and the choice h c = 12.53 obtained from the analytic study of the DPRM problem, 
which predicts a linear dependence exp(— lni?) = C[(h c /h) m — 1]. The circles are obtained with 
the projected cavity mapping, and a fitted value h c = 7.5. In both cases the exponent m is the one 
found from the DPRM analysis, m = 0.679. The straight line is a guide to the eye. 

the Bethe lattice, giving m = at the transition. Concerning the Monte-Carlo approach, let 
us point out in particular that the RSB effects happen at very low temperature, but they 
completely change the physics of the problem (for the problem with K = 3 and uniform 
couplings on [0, 1], we have found with the cavity- mean- field approximation that the RSB 
transition point on the ferromagnetic-paramagnetic transition line occurs at Tr ~ 0.087). So 
the numerical results obtained by extrapolation from a not-very-low temperature to T = 
might be questionable. 

VIII. SUMMARY 

We have proposed here a well denned mean-field scheme to study quantum ferromagnets 
in random transverse fields. This scheme consists in studying the problem on a random 
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regular graph with fixed connectivity K + 1. The problem at finite z is qualitatively distinct 
from the infinite z case, described by the naive mean field theory. It has a zero temperature 
quantum phase transition at a finite critical value of the width of the distribution of random 
transverse fields. 

The phase diagram has been studied by the RS quantum cavity method. A full study 
would involve a complicated population dynamics in terms of imaginary time spin trajec- 
tories, which would be extremely costly in terms of computing efforts. Instead we have 
resorted to a simple approximation, the cavity mean-field approximation, which parame- 
terizes the imaginary time spin trajectories by a single number, a local cavity longitudinal 
field. Within this approximation, the RS cavity method becomes a simple recursion rela- 
tion on these local fields. While this recursion looks structurally similar to the one found 
in classical spin systems, its physical content is very different. When one iterates L times 
the linearized recursion describing the phase transition, one obtains the partition function 
for a classical directed polymer in a random medium on a tree. This problem has a glass 
transition. The small disorder phase of the polymer describes the high temperature part 
of the ferromagnetic to paramagnetic phase transition, which is identical to what one gets 
in the naive mean field approach. The large disorder phase of the polymer describes the 
lower temperature part of the ferromagnetic to paramagnetic phase transition, which has 
a completely different behavior. This glass phase of the polymer, which displays replica 
symmetry breaking, actually describes the Griffiths phase of the quantum spin system in 
the neighborhood of its quantum phase transition. The results obtained with the cavity 
analysis fully agree with the known behaviour of the system in one dimension. On Bethe 
lattices with connectivity K + 1 > 3, we find a Griffiths phase, as in one dimension, but 
the critical point is not an infinite disorder fixed point: the exponent m which governs the 
tail of the local field distribution does not vanish at the transition, meaning that the critical 
exponent z does not diverge, it would be interesting to know if there is a critical dimension 
of finite dimensional problems beyond which such a behaviour appears. 

One should be aware of the fact that we have been using here a RS cavity method for 
the quantum spin system. Nevertheless, even within this RS method, the cavity equations 
map onto a classical problem which exhibits RSB. It would be very interesting, in particular 
when applying this formalism to spin glasses, to study the effect of RSB within the quantum 
cavity method itself. 
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